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Abstract 

Strong pionic decays of baryons are studied in a non-relativistic 
quark model framework via a convergent (p/E) expansion of the tran¬ 
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1 Introduction 


There is a current experimental and theoretical interest in the study of the 
baryon decays (see for instance [1]). Measurements which are planned at dif¬ 
ferent experimental facilities, TJNAF (CEBAF), MAMI, ELSA, GRAAL, aim 
to answer open questions about the baryon resonances, in particular having to 
do with the position of the Roper resonance in the baryon spectrum and with 
the absence of strong experimental evidence for many resonances around 1.8-2 
GeV. Both subjects are closely related to the study of the baryon decay ampli¬ 
tudes. Some authors have proposed that small couplings in the Nn production 
channel can provide an explanation for the missing states [2]. On the other 
hand the analysis of the corresponding decays is essential to any attempt to 
establish the true nature of the Roper resonance for which several alternatives 
have been proposed (hybrid three quark-|-one gluon state, first radial nucleon 
excitation...) Hence, the study of the Nn decay amplitudes becomes of great 
interest, as a test for the proposed models and in order to have a guide-line for 
future experiments. In this article we shall center our attention at the lowest 
energy decays for reasons that will become clear later on. 

In a previous paper [3] we put the emphasis in the construction of the 
effective transition operator for B B'n decays, making clear not only the 
very relevant known role played by the pion structure (otherwise said, by the 
qqq qq effective component of the baryon wave function) but also the need 
of taking into account relativistic effects. Here we go further in this analysis 
and show that ’better accounted’ relativistic corrections through a convergent 
{p/E) (instead of {p/m)) expansion does not spoil our results but instead 
conhrm in a more sound way our previous conclusions. 

For the sake of completeness we briefly review next the quark models used 
in ref. [3] to calculate the baryon wave functions which we shall employ for 
estimating baryon decay transition amplitudes. There are many other com¬ 
petitive models available in the literature. However for our purpose here, the 
comparison of a {p/E) expansion calculation against a {p/m) one, we do not 
feel our particular choice to suppose any serious limitation. 

2 The Quark Models. 

A precise description of the octet and decuplet baryonic spectrum (with and 
without strangeness) including the second excitation energy have been consis¬ 
tently obtained by solving the Schrodinger equation with a potential containing 
apart from the ’minimal ingredients’ (conhnement -f (coulomb -|- spin-spin) 
one gluon exchange interactions) a three quark phenomenological force [4]. 
Specihcally the expression used for the potential is: 
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where a, k, tq, Vq, mo are free parameters fixed from the spectrum and the the 
quark masses are chosen to get the baryon magnetic moments (see table 1). 
is the interquark distance and the a’s are the Pauli matrices. D is a constant 
to £x the absolute value of the nucleon mass to its experimental value. 

Such a model is adequate to solve altogether some endemic problems con¬ 
cerning the baryon spectrum, say, a unified description of the positive and 
negative parity states, the correct position of the Roper resonances (first ra¬ 
dial excitations) and the appearance of extra A negative parity states at a 
relatively low energy (A(l/2+) at 1900 MeV and A(5/2+) at 1930 MeV). 
Moreover, when combined with an improved transition operator, the B B'tt 
decay widths for the lowest energy non-strange resonances are in very good 
agreement (to the standard in the held) with the experimental data [3]. 

As a side effect, the splitting between the two lowest | and | states, is 
rather small for all the baryons studied (related to the absence of one gluon 
exchange tensor hyperhne or spin-orbit forces) and an unobserved proliferation 
of states would come out if we were to push the model further up the second 
excitation energy (maybe denoting the energy limit of the model description 
or maybe indicating the presence of missing states). 

We shall compare the results from the -|- P*-^^) potential model with 
the ones obtained with a two-body P^^^ potential model with parameters fitted 
to get an overall £t to the meson and baryon spectrum [5], the aim being to 
try to extract general features associated to ’two-body’ models and to make 
clear, if so, some bias of the results from our ’three-body’ approach. The small 
core size typical in spectroscopic models ((r^)^A ~ o.47 fm for P*^^)) makes 
the speed of the quarks be close to 1. As a matter of fact {p/m) > 1. In this 
sense a (p/E) expansion is more founded and incorporates automatically some 
relativistic corrections. 

Our results are expressed in terms of the width for each process. In all cases, 
the amplitudes have been calculated with the model wave functions, but we 
have used the physical masses for the kinematical factors. This is important 
for the potential model due to the discrepancies of its predictions with 
the experimental masses. 

Several mechanisms have been proposed to study the B B'tt decays. 
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3 Decay Mechanisms and Results. 

The elementary emission model considers the point-like pion emission by 
one of the qnarks of the baryon. The transition operator is obtained via 
the non-relativistic reduction of the qqTi interaction for which a pseudovec¬ 
tor form is used. The {p/m) order results have been published elsewhere [3]. 
The predicted widths of some low-lying non-strange baryons in the elemen¬ 
tary emission model are summarised in tables 2a and 2b, showing a general 
disagreement with experiment. Furthermore, the calculation of the quark ma¬ 
trix element h(p) 7 ^ 75 M(p) at the {p/m)‘^ order has been worked out [7] (we 
do not use the recalculated value of in order not to have an even bigger 
contribution), making clear the relevance of the corrections and the lack of 
convergence hence raising serious doubts about the adequacy to proceed to 
such an expansion. 

These problems can be partially palliated by proceeding to a {p/E) expan¬ 
sion {E = -|- m^) of the transition operator (for obvious reasons we shall 

work in momentum space). It is very instructive to do it order by order so 
that we can check the convergence of the {p/E) expansion and compare it to 
the {p/m) one. Thus the amplitude at {p/E) order, can be written as: 

= A)> = -n^S(PB - Pb' - k)A, (4) 
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where J,J' are the total baryon angular momenta, Pb,b' and k the three- 
momenta of the baryons and pion respectively, {m-,/) is the pion energy 
(mass), E (E’) refers to the energy of the incoming (outgoing) quark in the 
emission, T ($) stands for the momentum space (spin-isospin) wave function, 
(a*-^^) is the isospin (spin) Pauli matrix acting on the emitting quark, and 
fqqT, is the coupling constant in the qqn vertex and p^^ are the momenta 
associated to quark Jacobi coordinates. 
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The first term in the intermediate square bracket appearing in (5), arises 
(except for the ^ contribution) from the spatial components while all the 
other terms, proportional to arise from the time component of the current 
which is coupled to the derivative of the pion held. It is easy to verify that 
eq. (5) has well dehned properties (change in sign) when exchanging the role 
of the hnal and the initial baryon. Let’s also note that the pseudo-scalar Tiqq 
coupling leads to an expression similar to (5) where is replaced by the 
difference of the kinetic energies of quarks in the initial and hnal states, which 
is nothing else than the corresponding contribution to cUtt- 

Beyond the p/E order, typical terms like "^/{E + m), that appear in the 
Dirac spinors will be expanded as 


P 

E + m 


^(1 + ^ 

2E ^ 4^2 


+ ■■■), 


whereas a p/m expansion would give 
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As usual, we shall work in the center of mass system of the decaying baryon, 
i.e. Pb = 0. The only free parameter fqq^^ is hxed from the NNn form factor 
at zero momentum transfer {\k\ = 0). Additionally, the NNti form factor is 
predicted at \k\ 7 ^ 0 (both models do very well for low |A;|). fqqT, is determined 
to the order of the calculation from /Ar 7 V 7 r( 0 ). 

The results for the same set of widths are shown in tables 2a and 2b. By 
comparing columns four, hve and six (this one corresponding to the complete 
calculation for the pseudovector form) with column two and three, the better 
convergence of the {p/E) calculation becomes obvious for both quark models. 
In the p/E expansion, the {p/E)"^ correction make the result to tend smoother 
to the hnal one, while in the p/m expansion the {p/m)'^ correction either go 
in the wrong direction or make the result to overshoot the hnal one. Strictly 
speaking only for some cases: A(1232) —IYtt, iV(1520) —IYtt, iV(1535) — 
Avr, A(1650) —Avr and A(1700) —Att, a hrst order {p/E) treatment is 
justihed whereas for A(1440) —Ayr and A(1520) — Att it should be taken 
with caution. The situation becomes more dramatic when considering higher 
excitations for which we do not hnd a justihcation neither to a {p/m) nor to 
a {p/E) expansion. Let us note that for A(1440) —> Att, A(1600) — > Att 
decay amplitudes, the hrst order is suppressed by the orthogonality of the 
baryonic radial wave functions. The corrections (p/m)^ obviously break the 
orthogonality argument, hence the ehect is large for these transitions. Instead, 
for corrections (p/A)^, which are close to the unit operator, the orthogonality 
argument seems to work again, hence the smaller ehects. 
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Differences with the data require to look back at the reaction mechanism; 
in particular one can wonder whether the pion structure may play or not a 
signihcant role. 

The ^Pq quark pair creation model considers the creation of a qq pair 
in the hadronic medium that by later recombination gives rise to the outgoing 
pion. The effective transition operator is [ 8 , 9]: 


T = - 


^ / dpqdpq 

bl{Pq)d]{Pq) 


3jS{Pq + Pq) ^(110|m, -rn^y^ipq - Pq)2- j 


(9) 


where bl (dp are the i-quark (j-antiquark) creation operators, Pq {pq) is the 
three-momentum of the quark (antiquark) of the pair, is the solid harmonic 
polynomial, Z is the color-spin-isospin wave function of the pair, (110|m,-m) 
is the SU(2) Clebsch-Gordan coefficient and 7 is the coupling constant at the 
vacuum-gg vertex. Following our criterium given above, 7 is extracted from 

/vAr7r(0). 

Then the transition matrix element reads: 


(i?V|T|i?) =-37^(110|m,-m)/m (10) 
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where 
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where is the spin-isospin wave function of the pair and (<l>^) refers to 
the momentum space (spin-isospin) wave function of the pion. We choose for 
a gaussian form with the parameter R\ = 8 GeV“^ in order to reproduce 
the root mean square radius of the pion. 

The decay widths, calculated elsewhere [3], are compiled in tables 3a and 
3b. A comparison to the second column of tables 2a and 2b shows clearly the 
relevance of the pion structure. However, the predicted values differ very much 
from the experimental data. It is very illuminating to take the point-like pion 
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limit of the transition matrix element < B't:\T\B >. In ref. [3] this was done 
in conhguration space and compared to the {p/m) order expression snggesting 
a way to introdnce relativistic-like corrections in the scheme. The results 
obtained in this manner (we shall call it modihed ^Pq model (M^Pq)) shown in 
tables 3a and 3b, represent a rather amazing improvement of the predictions. 
Here we repeat the procedure in momentum space in order to get a ’more 
convergent’ {p/E) version of the modihed ^Pq model. 

In the point-like limit we get: 


(PV|T|P) 
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This expression reproduces (5) under the formal substitutions: 
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Eq. (13.a) points out an energy dependence (through the factor cny^) of the 
coupling constant 7 that can be associated to the boost from the rest frame of 
the pion to the decaying baryon rest frame. Moreover the expressions on the 
right hand side of (13.b) and (13.c) reduce to the left hand side in the very 
non-relativistic limit P ~ P' ~ m, (ci;^) 3 Pq ~ 2m. 

Hence we can go the other way around and use the right hand side of 
eq. (13) as the relativistic expressions of the left hand side, to modify the 
^Po effective transition operator. Technically the substitutions translate in 
the modihcation of the argument of the spherical harmonic plus an energy 
dependent factor: 


yT -\k 


3 1^6 + 


+ 



-^(1 + ^) + 



( 14 ) 
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Thus, the transition operator evidences symmetry properties expected from 
the elementary coupling of the derivative of the pion held to a current involving 
baryons. This ’relativized’ model that we shall call R^Pq from now on 
conhrm the good predictions the M^Pq version gave (the NNtt form factor is 
very similar to the one obtained with the usual ^Pq model). It is worthwhile 
to mention the excellent ht we get for the radial decays iV(1440) — Ntt, 
A(1600) — An, due entirely to the incorporation of the pion structure not 
having any need to consider any exotic nature for them. To this respect, it 
should be mentioned that the (p/P)^ and higher order corrections to them, 
although they were not small in the elementary emission model, are not very 
relevant compared to those arising from the pion structure. We get a good 
fit (exceptions are the iV(1535) —An and iV(1650) — An) even for some 
non-justihed hrst-order {p/E) treated cases maybe indicating the reabsortion 
of some relativistic effects through the only parameter. 

The 7V(1535) —Nn and iV(1650) —Nn decay widths are very badly 
described. For this several reasons can be pointed out. Higher order correc¬ 
tions are very signihcant as can be easily inferred from the tables. Besides, 
for these decays one should be aware of the presence of threshold effects [ 10 ] 
(the coupling for the Np channel is important). Additionally the lack of a 
tensor interaction at our model prevents mixing effects that could be relevant. 
To this respect we have evaluated the mixing that a tensor potential with 
the OGE parameter we have used would introduce. Though the mixing angle 
(■d ~ —40° is similar to the obtained with other models it does not translate 
in ours in a general improvement of all the A(1535) and A(1650) decays. 

In all cases, the effect of the emitted baryon width [11] and the hnal state 
interactions which have been neglected might play some role as well. Further¬ 
more an extended pion-quark vertex could be considered through a form factor 
multiplying (or 7 ). An estimation of its importance for a standard dipole 
type form factor [12] gives corrections to the widths of a 20 % at most. It 
is also worthwhile to realize the better predictions on the average that three- 
quark potential provide. They support the results for the baryonic spectrum 
and may indicate the convenience of the presence of high momentum compo¬ 
nents in the baryon wave function at low energies, though no strong conclusion 
should be derived from our restricted calculation. 

Other treatments of strong decays can be found in the literature within a 
^Pq scheme [13] getting a reasonable overall £t to a much wider range of data 
than considered here. However the adopted philosophy is quite a different 
one; we put the emphasis in the {p/E) expansion that allows the study of 
better relativistically controlled processes and show that a precise simultaneous 
description of the spectrum and the decays is feasible in some cases. 



4 Summary. 

We have performed a detailed analysis of strong pionic decay widths from a 
non-relativistic scheme, hrst following a {p/m) expansion approach, second 
using for the energy its relativistic expression {E = instead of the 

E K, m non-relativistic limit. A clear improvement of the convergence is ob¬ 
tained and a justihcation of a {p/E) treatment is obtained for some cases. The 
strong remaining discrepancies between our theoretical predictions and exper¬ 
iment has led us to analyse through a ^Pq pair creation model the influence 
of the qq content of the pion in the decay width and the need of introducing 
relativistic-like corrections into the effective transition operator. These rela¬ 
tivistic corrections have basically a kinematical origin An improvement of the 
{p/m) results for the examined cases is achieved in good agreement with exper¬ 
imental data even for some non-rapidly convergent processes. Let us mention 
in particular that the sizeable transitions between the N and A ground states 
and their radial excitations originate from the pion structure and not from the 
{p/m)‘^ terms in the operator as one may naively expect from the fact that the 
orthogonality of the states is of no relevance in this case. The validity of the 
argument is restored when the full series of higher order terms is considered. 
Being conscious of possible minor improvements that could be considered, as a 
general conclusion we might say that in order to do better for the pionic decay 
widths a more complete treatment of relativistic effects and the qq structure 
altogether in a consistent scheme seems to be unavoidable 
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Table captions 


Table 1 Fitted values of the parameters of the two-body and two+three 
body potentials. 

Table 2a Pionic decay widths and the fqqj^ coupling constant obtained with 
the elementary emission model (EEM) for the two-body quark potential 
model. Experimental data from [ 6 ]. 

Table 2b Pionic decay widths and the fqq.,^ coupling constant obtained with 
the elementary emission model (EEM) for the two+three body quark 
potential model. 

Table 3a Pionic decay widths and the 7 constant obtained with several 
versions of the model for the two-body quark potential model. 

Table 3b Pionic decay widths and the 7 constant obtained with several 
versions of the ^Pq model for the two+three body quark potential model. 
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Table 1 
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EEM EEM 
{p/m) (p/m)^ 


y(2) 

A(1232) ^ Nil 79.6 25.1 

7V(1440) ^ Att 3.4 177 

7V(1440) ^ Att 7.1 13.3 

A(1600) ^ Att 20.1 94.1 

A(1600) ^ Att 2.85 56.4 

A(1520)^A7r 61.8 21.5 

A(1520) ^ Att 78.0 95.3 

A(1535) ^ Att 240 494 

A(1535) ^ Att 9.7 4.25 

A(1650) ^ Att 47.9 135 

A(1650) ^ Att 12.4 5.49 

A(1700) ^ Att 4.07 1.22 

A(1700) ^ Att 383 612 

fqqiT 


0.602 0.602 


EEM 

(P/E) 

EEM 

ip/E)^ 

EEM 

All ip/E) orders 

Texp 

(MeV) 

74.3 

75.0 

75.8 

115-125 

0.007 

4.91 

13.1 

210-245 

5.4 

7.82 

9.30 

70-105 

31.0 

15.4 

7.98 

35-88 

0.72 

5.08 

9.30 

140-245 

62.3 

60.6 

57.7 

60-72 

24.3 

39.2 

45.1 

18-30 

24.8 

76.5 

101 

53-83 

9.93 

10.4 

10.3 

<1.5 

2.49 

14.5 

20.9 

90-120 

12.6 

13.0 

12.8 

4-11 

3.65 

3.43 

3.24 

5-15 

112 

190 

222 

81-393 

0.619 

0.743 

0.782 
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y(2) ^ y(3) 


EEM EEM 
{p/m) (p/m)^ 


A(1232) ^ Nil 72.1 4.18 

7V(1440) ^ Nil 0.17 690 

7V(1440) ^ Att 17.6 43.1 

A(1600) ^ Att 94.1 226 

A(1600) ^ Att 0.10 70.4 

A(1520) ^ Att 22.3 15.3 

A(1520) ^ Att 56.1 67.4 

A(1535) ^ Att 149 259 

A(1535) ^ Att 8.3 5.78 

A(1650) ^ Att 24.9 60.3 

A(1650) ^ Att 9.95 6.35 

A(1700) ^ Att 1.43 0.75 

A(1700) ^ Att 220 347 

fqqiT 


0.604 0.604 


EEM 

(P/E) 

EEM 

ip/E)^ 

EEM 

All ip/E) orders 

Texp 

(MeV) 

67.4 

71.7 

75.7 

115-125 

0.04 

4.01 

15.3 

210-245 

16.7 

24.5 

30.9 

70-105 

84.9 

64.1 

50.4 

35-88 

0.006 

2.08 

5.89 

140-245 

25.9 

29.3 

31.5 

60-72 

20.4 

36.2 

44.8 

18-30 

20.7 

60.7 

82.8 

53-83 

9.25 

11.3 

12.6 

<1.5 

2.54 

11.0 

15.8 

90-120 

10.4 

12.3 

13.6 

4-11 

1.29 

1.41 

1.52 

5-15 

79.7 

148 

185 

81-393 

0.624 

0.780 

0.852 
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y(2) 3p^ M3Po R^Po Te^p (MeV) 


A(1232) ^ Nil 167 

7V(1440) ^ Nil 452 

A(1440) ^ Att 66.5 

A(1600) ^ Nn 19.8 

A(1600) ^ Att 255 

A(1520) ^ Nn 268 

A(1520) ^ Att 532 

A(1535) ^ Att 429 

A(1535) ^ Att 28.1 

A(1650) ^ Att 49.1 

A(1650) ^ Att 49.3 

A(1700) ^ Att 11.5 

A(1700) ^ Att 1643 

—i7 7.02 


88.6 

83.9 

115-125 

114 

73.5 

210-245 

27.6 

20.7 

70-105 

2.1 

0.27 

35-88 

62.0 

41.9 

140-245 

95.1 

92.0 

60-72 

45.9 

17.1 

18-30 

49.2 

0.18 

53-83 

15.3 

15.7 

<1.5 

6.14 

1.13 

90-120 

20.9 

20.9 

4-11 

3.49 

3.15 

5-15 

228 

114 

81-393 


8.20 8.29 


Table 3a 
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y(2) + y(3) 3p^ M3Po R^Po Te^p (MeV) 


A(1232) ^ Nil 210 

7V(1440) ^ Nil 1076 

7V(1440) ^ Att 228 

A(1600) ^ Nn 0.53 

A(1600) ^ Att 498 

A(1520) ^ Nn 319 

A(1520) ^ Att 999 

A(1535) ^ Att 464 

A(1535) ^ Att 74.0 

A(1650) ^ Att 44.2 

A(1650) ^ Att 109 

A(1700) ^ Att 11.2 

A(1700) ^ Att 2417 

—i7 9.78 


112 

106 

115-125 

307 

236 

210-245 

116 

106 

70-105 

2.5 

8.82 

35-88 

121 

93.5 

140-245 

105 

100 

60-72 

75.2 

34.4 

18-30 

44.1 

0.52 

53-83 

36.3 

37.7 

<1.5 

4.01 

0.77 

90-120 

45.7 

43.1 

4-11 

3.25 

2.83 

5-15 

323 

217 

81-393 


11.4 11.6 


Table 3b 
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